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t/3 ' Abstract 

Given a collection {Zi, Z2, ■ ■ ■ , Zm} of n-sided polygons in the plane 

and a query polygon W we give algorithms to find all Ze such that W = 

f{Zi) with / an unknown similarity transformation in time independent 

^ ■ of the size of the collection. If / is a known affine transformation, we show 

\l ' how to find all Ze such that W = f{Ze) in 0{n + log(rn,)) time. 

O^ , For a pair W, W' of polygons we can find all the pairs Zi, Zgi such that 

'**'/ ■ W — f{Ze) and W' = f{Zii) for an unknown affine transformation / in 

0{m + n) time. 

^^ , For the case of triangles we also give bounds for the problem of match- 

^3 ' ing triangles with variable vertices, which is equivalent to affine matching 

Cn ' triangles in noisy conditions. 



1 Introduction 

The problem of matching point sets under similarities and affine transforma- 
tions has been addressed by matching sets of triangles, extracted in a consisting 
way, from the point sets [S]. In this work we propose a technique to matching 
arbitrary polygons under affine transforms, and as a particular case triangles in 
noisy conditions. 

We will identify points {x, y) in the plane with the corresponding complex 
numbers z = x + iy. A polygon in the plane will be an ordered set of points, or 
complex numbers. Since the set is arbitrary, self intersection and crossings are 
allowed. 

An afhnc transformation / : R^ — > M^ can be (uniquely) written as 

f{z) ~ az + j3z + "f 
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where a,/3,j € C and jap — |/?p ~ det/ y^ 0. Here z stands for the complex 
conjugated of z. 

Given polygons Z ~ (zi, 2:2, . . . , Zn) and W — {wi,W2, • • • , w„) with z^ and 
We complex numbers, the problem consist in determining if there exists an affine 
transformation / such that Z = f{W). Since the affine transformation have 
three complex parameters, it is enough to find two corresponding triples of 
consecutive points in both polygons. A naive procedure will be to fix a triple 
in Z and try all the shifts in W to find the correspondence. This takes 0{n) 
operations. 

Now assume we have a given collection of polygons Zi, Z2, . . . , Zm and a 
query polygon W, and want to know which of the Zg arc affine images of W . 
Using a sequential approach and the simple procedure above the solution can 
be found in 0{mn) operations. In general, without an index, the complexity 
will depend linearly on the number of polygons in the collection; multiplied 
by the complexity of an individual match. We will show how to improve this 
complexity using invariants. 

2 Invariants 

We split the invariant in two cases because similarity transformations (rotations, 
translations, scaling and compositions between them) have relevance per se. 

2.1 A Similarity Invariant for Polygons 

Let j and n be integers, with n > 3, and let p : {1, 2, . . . , n} — > {1, 2, . . . , n} be 
a permutation. Consider functions (pn.p, fnj '■ C" — > C U {00} given by 



ipp{zi,...,Z.a) ^ 
(^„j(zi,...,Z„) 






En X —JK^ 
fe=l ■^ri Zk 



where A„ = e^'^*/". 

The function ipp is well-defined except on 

n n 

AA=.{(zi,...,z„)eC" : ^APWzfe = = ^A-P«z,}, 

fe=l k=l 

but A/" is a n — 2 dimensional vector subspace with null measure in C" . Similarly, 
ipn,j is well-defined except on a n — 2 dimensional vector subspace. 

Notice that ip-p is invariant under the action of similarity transformations on 



polygons with n vertices: 



ipp {azi + 13, az2 + /3, . . . , az„ + /?) 






(^p(zi,Z2,---,Zn)- 



It is also clear that Lpn^j is invariant under the action of similarity transfor- 
mations: 

LPrL,]{aZi + P,aZ2 + l3,...,aZn +/?) = ipn.jizi,Z2,...,Zn). 

2.2 An AfRne Invariant for Polygons 

In this section, (p denotes either cpn.p or (pn.j- Let ( — (p{zi, Z2, ■ ■ ■ , Zn) and ^ = 
'/'(/(-^i); /(-^z), ■ ■ ■ , f{zn)), where / is an affine transformation. By making x = 
I]fc=i ^n Zk and y = J2'k=i A"*" -2;^ ii if = ipp, or by making a; = Y.k=i ^n Zk 
and y = X]fe=i ^n''' ^k if V? = '/'nj, we obtain 
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(2.1) 



That is, the number = — does not depend neither on the n-agon (zi, . . . , 2„) 

nor on the numbers n and j, nor on the permutation p; it just depends on /. 

Please notice that similarity invariance in section 12.11 can be obtained from 
the above taking (5 — 0. 

Remark 2.1. The numerator and denominator involved in the definition of ipnj 
are coefficients of the polygon (zi, . . . , z„) appearing when it has expressed in 
certain basis of C" , namely the basis of star-shaped polygons (see [2] , [H proof 
of Proposition 3]). 

Remark 2.2. Equality ()2.ip is analogous to a well-known result from Teichmiiller 
theory: if we deform the complex structure [T] of a torus by an affine transfor- 
mation /, then the Teichmiiller distance between [T] and [/T] does not depend 
of r, it just depends on / (see [31 chapter V.6] for instance). 

2.3 Shifts 

A polygon can be enumerated in n different cyclic ways or shifts. Notice that 
V3„j(zi, Z2, ■■■TZn) = \'^(pn.j{z2, Z3, . . . , z„ , zi) , heuce we have the equality of 
potencies 



Vn,j(zi,Z2, . . . , Z„)" — ipnj{zi, Z|+i, . . . , Z£+„_i , Zf+n)*" 

for any £, where the subscripts are taken mod n. 



(2.2) 



Remark 2.3. In what follows we use functions of the form (pn,j because (pf, does 
not satisfies an identity of the form (|2.2p for a general permutation p. Each 
one of the instances of (pn.j, where j ranges in {1, 2 . . . , n — 1}, can be used in 
the below algorithms for matching polygons. Since they will be used as hashing 
functions, collisions can be avoided by using multiple (p„j instances. 

2.4 An Index for Matching Polygons 

Assume that a collection of different polygons Zi, Z2, . . . , Z,„ of n edges is given. 
By a preprocessing step we compute pairs {i, ipnj{Zi)"). Assume that a query 
polygon W is given and that the objective is to find all the polygons in the 
collection such that W — f{Zi) for some unknown similarity transformation 
/. Using the result from section [^751 this corresponds to all the polygons such 
that (pnj{Zg)" = ipn.j{W)^, assume there are R of them. They can be found in 
0{n + R) operations. Please notice that this search can be implemented using 
hashing in time independent of the size of the collection. 

If the parameters of the affine transformation are known, then all the match- 
ing polygons such that W = f{Zi) can be easily found. Using the result 
from section 12.21 we compute C, = (pnj{W). If the polygons match, then 



IC.-CI 



/3 



, where ^e = fn.jiZe)- We can retrieve candidate objects in 

sublinear time using a spatial access method (e.g. kd-trees), or a metric in- 
dex [1]. The corresponding R candidate objects in the complex plane will be 
arranged as a circle surrounding the query polygon W. Using an index, e.g. 
kd-trees, the complexity bound of the search will be 0{n + log(TO) -|- Rt) with 
n the number of edges of the polygon, and m the number of polygons in the 
collection and r the cost of checking the matching between a pair of polygons. 

2.5 An Index for Matching Pairs of Polygons 

If the affine transformation is unknown, the search space cannot be bounded 

using invariant 12.21 However, the problem can be solved if we have a pair of 

polygons. Let W and W be two n-sided query polygons, we want to obtain 

all polygons Zg, Zgi such that W — f{Zi) and W — f{Zg') for some unknown 

affine transformation /. We compute C, — (pn,j{W) and C = (pn.j{W). Let 

\F — C\ \f — C'\ 
" — and ry^ — = — with ^^ — Lpnj{Zi). Abusing the notation. 



let {Zi} = {ije} and {Z'^} — {r/^}. Using the result from section [521 what we 
need is the intersection of the two sets above {rjg} n {rf(}. The polygons in the 
intersection are candidates for matching. The algorithm above needs 0{n + m) 
operations. To finish the procedure, 0{R) match verifications of the R polygons 
in the intersection are needed. 



2.6 Collisions 

For two unrelated polygons Wi,Wi it is possible that they collide, i.e. ipn,j{Wi)^ - 
iPnj{Wi)" without being an affine related. This increase the R candidate objects 
to be reviewed in algorithms above. To decrease the number of candidates we 
can use multiple instances of ipuj- Let (pn,i, <Pn,2, • • • , ¥'«,«-! be the collection of 
instances of the affine invariant functions, and Li,L2, . . . ,in-i the respective 
lists of candidates. The true matches should be in the intersection. Notice that 
it is not necessary to use all the (pn.j functions, only a subset of them. 

3 Noisy Polygons 

A slightly more general setup is when there is an unknown noise function in the 
matching. The image of the query polygon is an affine transformation plus noise, 
namely W — {f{zi + Azi), f{z2 + AZ2), . . . , f{zn + Azn))- We have bounded 
the difference for the case of triangles, this is explained below. We were not able 
to give a tight bound for general polygons. To simplify notation, A3 = e^'^*/'^ 
and (ps^i will be denoted respectively by A and (p in what follows. 

3.1 Triangles with Variable Vertex 

For the similarity class of a triangle (21, ^2, Z3) we can choose a representative 
of the form (0,1, r). The complex number t = t(zi, 22,23) is determined as 
T = (z3 - zi)/{z2 - zi). Let us define M{t) == (A^ + r)/(A + r). This is the 
unique conformal biycction M : C U {00} — > C U {00} such that 

ip{zi,Z2, Z3) = M(t(zi, Z2, Z3)). (3.1) 

The order 3 'rotation' R{t) = 1/(1 — r) of CU {00}, with fixed point — A^, cor- 
responds to cyclic relabelling of the vertices of a triangle, because the triangles 
(0,1, r), (-R(t), 0,1) and (l,i?^(T),0) arc similar. We use the term 'rotation' 
because M o Ro M^^(^) = A^ is a rotation by angle 2tt/3. 

We shall proof an estimation for perturbations of an equilateral triangle. For 
a real number < r < \/3/6 we consider the polydisc 

C/r = {(21,^2,^3) eC^ : |zi-A| <r,\z2-X^\ < r, |z3 - 1| < r}. (3.2) 

Triangles (zi, Z2, z^) e Ur can be seen as perturbations of the equilateral triangle 
(A,A2,1). 

Lemma 3.1. 

9 - 4^/3r + I2r^ - (3 + 2VSr)\/9 - 20V3r + I2r^ 



\ip{zi,Z2,Z3)\ < 



\ 9 - Ay/Sr + 12r2 + (3 + 2V3r)V9 - 20V3r + 12r2 



-^V3r+|r2 + 0(. 



for any (21,^2,2:3) e Ur 
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Figure 1: The Apollonian circles Cki and Ck^, where Ki = (\/3 — 2r)/(\/3 + 2r) 
and K2 = 1/Ki, which have respectively diameter a;iX2 and xj,X4^ with xi = 
2r - v^ \/3 - 2r \/l + 2r VS + 2r 
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Here 0{r^) stands for the product of r^ and a convergent power series in r. 

Proof. First we claim that if we see triangles (zi, Z2, 2:3) € Ur as triangles in the 
form (0, 1, r), then r lies in the non-shaded region showed in fig. [21 which will 
be described below. To proof the claim, note that the length of the longest edge 
of a triangle (zi, Z2, 23) with vertices in Ur is at most \/3 + 2r, and the length 
of the shortest edge is at least -\/3 — 2r. It follows that 



A/3^2r 

%/3 + 2r "" \zj~ zi\ " V3-2r 



< 



Zfel 73 + 2r 



(3.3) 



with {j, fc, 1} = {1, 2, 3}. ApoUonius' theorem says that if iC > 0, then the locus 
Ck = {t G H : |t — 1| = i^|T|} is an euclidean circle centered at the real axis. 
Since (A, A^, 1) is a positive oriented triangle, we have that the imaginary part 



of T is positive. Inequalities (|3.3p imply that r can not lay in the shaded regions 
of fig. [U and T can not lay in the rotations by R and i?^ of the shaded regions. 
Hence r should lay in a curvilinear hexagon 'centered' at — A^ (see fig. [2]). The 
claim is proved. 

To conclude the proof of the lemma, consider the vertex 



^ 1 . \/9 - 20\/3r + 12r2 

C = \- I 7= 

2 2(V3 + 2r) 

showed in fig. [S] The points i?(C), R^{0 are also vertices of the curvilinear 
hexagon. The transformationation M maps the circle through C,, R{C)j -^^(C) 
onto a circle of radius 



|M(C)| 



centered at the origin. 

We have hm^^^g^g- ( = 1/2 e 
triangle. 



9 - 4v/3r + 12r2 - (3 + 2V^r)\/9 - 20\/3r + 12r2 
^ 9 - 4V3r + 12r2 + (3 + 2V3r)\/9 - 20A/3r + 12r2 

and (0, l,r) with t G M is a degenerated 
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Figure 2: Curvilinear hexagon centered at — A'^ with a symmetry group 



The first two coefficients of the Taylor series result from a straightforward 
calculation. D 

Now we deal with perturbations of a general triangle {zi,Z2,zz). 
Notice that the unique afSne transformation /(z) = az -\- /3z + 7 that takes 
the equilateral triangle (A, A^, 1) to a non-degenerated triangle (zi, Z2, z^) is 

X'^ZI+XZ2 + Z3 Xzi + X'^ Z2 + Z3 _ Z1+Z2 + Z3 

3 " + 3 "+ 3 • (^-'^ 

In order to give a geometric interpretation to the affine transformation f{z) = 
az -|-/3z + 7, we set a — jaje^'^^ /3 — \j3\e^'^'^ and z — re*^. Then az + /3z + 7 is 
equal to 



re 



i(0l+02)/2 



(l"l + l/?l)' 



I sm 



- '/'2 



e 



-1- 



Therefore / maps a circle of radius r onto an ellipse with major axis of length 
2r(|Q;| + |/3|) inclined at angle {(pi + ^2)/2 with the real axis, and minor axis of 
length 2r||Q;| — |/3||. It is convenient introduce Ez{pi,p2,0) to denotes the open 
(i.e. without boundary) ellipse with major axis of length pi inclined at angle 
9 with the real axis, minor axis of length p2, and centered at the point z. The 
proof of the next result is immediate from the above geometric interpretation 
for affine transformations. 

Lemma 3.2. Given < r < \/3/6 and a non- degenerated triangle, we set 



A^zi + Az2 + Z3 



be">'^ = 



Azi + A^Z2 + Z3 



with a,b > 0. Let V = {(wi, 102,^3) G C^ : Wj € Ezj{pi,p2,0) for j = 1,2,3} 
where pi — 2r(a + b), p2 — 2r\a ^ b\, 6 = {(pi + 02)/2. Then the affine trans- 
formation such that takes (A, A^, 1) to (zi, Z2, 23) maps triangles with vertices in 
Ur (see eq. (|3.2[) ) onto triangles with vertices in V . 

By the proof of lemma lXTl triangles {zi, Z2, 2:3) G Ur are similarity equivalent 
to triangles (0, 1, t) with r in a region inside the circle C through (, R{0' ^^(C)j 
where 

1 . V9 - 20V3r + 12r2 

2^' 2(x/3 + 2r) ' 

and this circle is mapped by M to a circle centered at the origin. Circles centered 
at the origin are preserved under multiplication by —1, therefore C is preserved 
under Af~^(— Af(T)) == (r — 2)/(2t — 1). Hence C is a circle with diameter 
between C and 



2 2V9-20\/3r + 12r2 

The hypothesis in the next result are the same as in the previous lemma. 

Theorem 3.3. Given < r < \/3/6 and a non-degenerated positively oriented 
triangle (zi, Z2, Z3), we set 

,4,, ^ A^zi + Az2 + 2^3 , .0, ^ Azi + A^Z2 + Z3 

3 ' 3 ' 

with a,b > 0. Let V = {(lUi, u'2, 1^3) G C^ : Wj G E)zj{pi, P2^Q) for j = 1,2,3} 
where pi = 2r(a + 6), p2 ~ 2r\a — 6|, 9 — (0i + (/)2)/2. T/ie?i i/ie triangles with 
vertices in V are similar to triangles (0,1, r) with t in a region inside of the 
ellipse 

^ (2V3p{a + b) 2V3p\a^b\ 0i + 02 , A 

Ew —: 7—, —: 7-, 7y arg(Z2 - ^1 ) 

\_ \Z2-Zi\ \Z2-Zi\ 2 J 

where 

1 A/3(zi+Z2-2z3)(9 + 12r2-4V3r) 



w ^ — — 



2 6(z2 - zi)(V^ + 2r) ^9 - 20V3r + 12r2 

8V3r 

p — I. 

(V3 + 2r) V9 - 20%/3r + 12r2 

Proof. We saw above that triangles (zi, 22,23) € Ur are similar to triangles 
(0, 1, r) with T inside the circle C of radius 

le - CI SySr 



{VS + 2r)\/9 - 20V3r + 12r2 

centered at (C + ^)/2; but now we are concerned with triangles with vertices in 
V. 



Let / be an affine transformation of the form f{z) — az + fJz + j on t. 
We have (/(0),/(1),/(t)) = (7, a + /3 + 7,ar + /3r + 7), which is similar to 
(0, 1, (ar + fiT)/{a + /3)). That is, the afhne transformation /, acting on the 
three vertices of the triangles, translates to the affine transformation /i(t) = 
(ar + /3T)/(a + /?) acting on the vertex r of (0, 1, r). We apply this remark to 
the transformation / = (72 o <7i where gi and 52 are the affine transformations 
which maps respectively (0, 1, — A^) to (A, A^, 1) and (A, A^, 1) to (2:1,22, 23)- For 
this f{z) = (A - l)ae*'^i z + (A^ - l)&e*'^^ 2 + ae*"^! + 6e* + (21 + 22 + 23)73, the 
corresponding /i(t) is 

hir) = T + -^ T. 

22 - 21 22 - 21 

Finally, the circle C is mapped by h to the ellipse 

^ (2V3pia + b) 2V^p\a-b\ 0i + ,^2 , A 

\^ I22-21I |22-2i| 2 J 

,fC + ^\ 1 ^/3(2l+22-223)(9 + 12r2-4^/3r) 

where w = a — — — = . □ 

V 2 y 2 6(22-zi)(v^+2r)V9-20\/3r + 12r2 

Remark 3.4. To obtain a bound for the image under ip of triangles with vertices 
in V, it is sufficient apply the transformation M (see equation p.ip ) to the 
ellipse 

^ /2V3p(a + 6) 2y3p|a-5| (j>i + (j>2 , \ 

y I22-21I I22-21I 2 J 

of theorem 13.31 
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